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ABSTRACT 



We have used the Australia Telescope National Facility Mopra 22m millimetre 
telescope to search for emission from the 6_2-7_i E (85.5-GHz) and 72-63 A~ (86.6- 
GHz) transitions of methanol. The search was targeted towards 22 star formation 
regions which exhibit maser emission in the 107.0-GHz 3i-4o A + methanol transition, 
as well as in the 6.6-GHz 5i-6o A + transition characteristic of class II methanol maser 
sources. A total of 22 regions were searched at 85.5 GHz resulting in 5 detections, of 
which 1 appears to be a newly discovered maser. For the 86.6-GHz transition observa- 
tions were made of 18 regions which yielded 2 detections, but no new maser sources. 
This search demonstrates that emission from the 6_a-7_i E and 72-63 A~ transitions 
is rare. Detection of maser emission from either of these transitions therefore indicates 
the presence of special conditions, different from those in the majority of methanol 
maser sources. We have observed temporal variability in the 86.6-GHz emission to- 
wards 345.010+1.792, which along with the very narrow line width, confirms that the 
emission is a maser in this source. 

We have combined our current observations with published data for the 6.6-, 12. 1-, 
85.5-, 86.6-, 107. 0-, 108.8- and 156.6-GHz transitions for comparison with the maser 
model of Sobolev & Deguchi (1994). Both detections and nondetections are useful 
for setting limits on the physical conditions in star forming regions which contain 
methanol maser emission. This has allowed us to estimate the likely ranges of dust 
temperature, gas density, and methanol column density, both for typical methanol 
maser sources and for those sources which also show 107.0-GHz emission. The gas 
temperature can also be estimated for those sources exhibiting masers at 85.5 and/or 
86.6 GHz. 
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1 INTRODUCTION 

More than 30 different transitions of methanol have been 
observed to exhibit interstellar maser emission. The maser 
transitions are empirically classified into two groups (named 
class I and class II) on the basis of their association with 
high- mass star formation regions IMenten 1991at . Class II 
transitions are closely associated with signposts of high- mass 
star formation such as ultra-compact Hll regions, OH masers 
and strong far infra-red emission. Class I masers occur off- 
set from star formation regions and are believed to be colli- 
sionally pumped at the interface between outflows and the 
parent molecular cloud IMenten 1991all . 

In contrast to OH and H2O masers, class II 



methanol maser emission is only found associated 
with high-mass star formation, making it a useful 
signpost of such regions. In particular the 6.6-GHz 
5i-6o A + transition has been detected towards more 
than 550 sites throughout our Galaxy (Mentcn 1991b 
|MacLeod, Gaylard fc Nicolson 1992| ICaswell et al. 1995al 
|Ellingsen et al. 1996| IWalsh et al. 19971 |Slysh et al. 19*9*5] 
|Szymczak et al. 2002^ . More than 120 of the regions 
which exhibit 6.6-GHz maser emission also exhibit 
maser emission in the 12.1-GHz 2o-3_i E transition 
JCaswell et al. 1995bl . Approximately 175 known 6.6- 
GHz maser sources have been searched in the 3i-4o A + 
transition at 107.0 GHz, with maser emission detected 
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towards 25 regions (TVal'tts et al. 19951 IVaFtts et al. 19991 
ICaswell et al. 20001 IMinier fc Booth 200211 . In contrast, the 
scale of searches and the number of detections for the other 
20 or so class II masing transitions have been very limited. 

Sobolev & Deguchi ( 1994;) developed the first theoret- 
ical model to successfully reproduce the high brightness 
temperatures that characterise strong class II methanol 
masers (referred to hereafter as the SD model). In the SD 
model the masers are pumped by infra-red radiation from 
warm dust close to the star forming region, where the gas 
phase abundance of methanol is greatly enriched following 
grain mantle evaporation. The SD model finds that the 6.6- 
and 12.1-GHz transitions produce strong maser emission 
over a wide range of physical conditions, providing an 
explanation as to why maser emission is common in these 
two transitions jSobolev, Cragg fc God frey 1997liJ. It also 
predicts weaker maser emission from a variety of methanol 
transitions ( |Sobolev, Cragg fc Godfre y J^97b J, including 
all the transitions discovered previously IWilson et al. 19841 
IWilson et al. 19851 |Haschick, Baan fc Menten 19891 

IWilson et al. 19931 [ Slysh, Kalenskii fc Val'ttsl995T and 
subsequently (IVaTtts et al. 19991 |Cragg et al. 2001^ . 

There is considerable evidence that many of the 
class II methanol maser transitions are spatially co- 
incident, arising from the same region. The 6.6- and 
12.1-GHz methanol masers have been shown to be co- 
incident at the milliarcsecond level towards a num- 
ber of sources (IMenten et al. 19921 IJNorris et al. 19931 
|Minier, Booth fc Conway 2000[ l . The higher frequency weak 
maser transitions typically align in velocity with 6.6- and 
12.1-GHz spectral features, but definitively proving that 
they are coincident is more difficult. However, recent BIMA 
observations of the 86.6- and 86.9-GHz masers in W3(OH) 
show that they are coincident with the strongest 107.0 
GHz masers to within the O'.'l accuracy of the observations 
(Sutt orTet al. 200Tt . These observations support the propo- 
sition that the class II methanol maser transitions are ex- 
cited simultaneously, which if true means that their obser- 
vation provides tight constraints for maser pumping models. 

Many transitions become masers simultaneously in the 
SD model, so that when masers at several frequencies co- 
incide, the best fitting model conditions can be estimated 
via multi-transition analysis. This provides a probe of the 
physical conditions in the maser region. The maser transi- 
tions which are strongly inverted only over a small range of 
model conditions are potentially much more useful for set- 
ting constraints on physical conditions than are the stronger, 
more ubiquitous maser transitions. For example in the mod- 
eling of Sobolev et al. the 85.5-GHz transition is predicted 
to be strongest for gas temperatures less than 50K and gas 
densities in the range 10 6 - 10 s cm , while the 86.6-GHz 
transition favours higher temperatures and lower densities. 
The observations of W3(OH) by Sutton et al. (1217011 de- 
tected maser emission from the same location in the 86.6-, 
86.9- and 107.0-GHz transitions, but only thermal or quasi- 
thermal emission from the 85.5-GHz transition, consistent 
with a warm gas model. In contrast, towards the southern 
star forming region 345.010+1.792, 85.5- 86.6-, 86.9-, 107.0- 
and 108.1-GHz masers have all been detected from the same 
velocity range (IVal'tts et al. 19991 |Cragg et al. 20"oT) . 

Previous searches for the 85.5- and 86.6-GHz methanol 
maser transitions have been made towards 30 and 35 



sources respectively i Cragg et al. 200l| ISutton et al. 20011 
IMinier fe"F jooth 2002 . In this work we report the results of 
observations of a further 15 and 10 sources not previously 
searched for 85.5- and 86.6-GHz methanol maser emission 
respectively. Both weak maser and quasi-thermal emission 
have been detected, and it is not always easy to distinguish 
between the two in single dish observations. The 85.5-GHz 
emission is thought to be of maser origin in 345.010+1.792, 
9.621+0.196, 29.95-0.02 and DR21(OH), while the 86.6- 
GHz emission shows maser characteristics towards W3(OH), 
345.010+1.792 and W51-IRS1. The search described here 
forms part of a program to observe a number of the rarer 
methanol maser transitions in a consistent sample of sources, 
selected on the basis of their 107.0-GHz maser emission. Of 
the 25 known 107.0-GHz masers, 22 can be observed with the 
Mopra 22-m telescope, and we have searched for 85.5- and 
86.6-GHz methanol masers towards those 107.0-GHz masers 
not targeted in previous surveys. Since only 4 of the 7 pre- 
vious detections were associated with 107.0-GHz emission, 
this strategy will provide useful information for modelling 
the target sources, but does not encompass all candidates 
for new 85.5- and 86.6-GHz emission. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were made between 2002 July 31 and 
August 13 using the Australia Telescope National Facility 
(ATNF) 22-m millimetre antenna at Mopra. At 86 GHz the 
FWHM antenna beamwidth was 39". The sources were ob- 
served in position switching mode, with 229 seconds spent 
at the onsource position and 229 seconds offset by -1' in 
declination. For the majority of sources this procedure was 
repeated 10 times giving a total onsource integration time 
of 38 minutes, which typically yielded an RMS of 0.7 Jy. 
Two different transitions were observed, the 6_2-7_i E and 
72-63 A~ transitions for which the adopted rest frequen- 
cies were 85.568084 GH z and 86.615578 GHz respectively 
JTsunekawa et al" 199511 . The data were collected using a 
2-bit digital autocorrelation spectrometer configured with 
1024 channels spanning a 64-MHz bandwidth. For an observ- 
ing frequency of 86 GHz this configuration yields a natural 
weighting velocity resolution of 0.26 km s _ , or 0.44 km s _1 
after Hanning smoothing. 

The antenna pointing was checked approximately every 
two hours through observations of nearby strong SiO masers 
and showed an RMS scatter of the order of 10". The sys- 
tem temperature was determined by inserting an ambient 
temperature load which was assumed to have a tempera- 
ture of 295K. The measured values varied between 210 and 
300K during the observations depending upon the weather 
conditions and elevation. This method of calibration also 
corrects for atmospheric absorption (Kutncr & Ulich T98ll 
and taking into account pointing inaccuracies the absolute 
flux density scale of the observations is believed to be ac- 
curate to 20 per cent. The amplitude scale was converted 
from antenna temperature to Jy assuming a sensitivity of 
30 JyKT 1 . The data were processed using the SPC reduction 
package. Quotient spectra were formed for each on/off pair 
of observations, which were then averaged together, a poly- 
nomial baseline fitted and subtracted and the velocity and 
amplitude scale calibrated. 
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An erroneous observing procedure meant that at each 
of 85.5- and 86.6-GHz three sources were observed with the 
telescope offset by +20" in right ascension (approximately 
half the FWHM of the antenna beam at this frequency) from 
the nominal position. The upper limit listed in Tables 0& 
121 for the sources observed at an offset position has been set 
to 9 (3 x 3) times the RMS noise level in the final spectrum 
as a conservative estimate of the detection limit. 

3 RESULTS 

3.1 6-2-7_i E (85.5 GHz) 

A total of 22 sources were observed at 85.5 GHz and these 
are listed in Tabled With the exception of Orion KL each 
of these sources shows maser emission in the 6.6- and 107.0- 
GHz transitions. The observations of Orion KL were used to 
test the system performance and to aid comparison of the 
calibration with the results of Cragg et al. <!2001t and Minier 
& Booth Emission from the 85.5-GHz 6_ 2 -7-i E 

transition of methanol was detected towards 5 sources. The 
spectra of these sources are shown in Fig. Three of the 
line profiles are broad (Orion KL, NGC6334F and 23.440- 
0.182) and so the emission is thermal, or quasi-thermal in 
origin. The emission from 345.010+1.792 is very strong and 
narrow, and as it aligns with the maser emission at other 
frequencies and is offset from the velocity range of the ther- 
mal emission in this source, it is clearly also masing in this 
transition. The 85.5-GHz emission towards 328.808+0.633 is 
narrow, but much weaker than 345.010+1.792 and so from 
these observations we cannot be certain that it is a maser in 
this transition. However, the velocity of the emission coin- 
cidences exactly with the velocity of the 107-GHz methanol 
maser in this source, adding support to the argument that 
328.808+0.633 is an 85.5-GHz maser. The observation of 
9.621+0.196 with 3-cr=1.5 Jy was not sufficiently sensitive 
to confirm the previous detection of this source with a peak 
flux density of 1.2 Jy (Cragg et al. 2001 ). For modelling we 
have used the flux density measured by Cragg et al. rather 
than the current observations. 

Gaussian profiles have been fitted to each of the sources 
where emission was detected and parameters are listed in 
Table [I] 



3.2 7 2 -6 3 A" (86.6 GHz) 

The 86.6-GHz 72-63 A - transition was detected towards 
3 of the 18 sites observed (see Table The spectra of 
sources where emission was detected (with the exception 
of NGC6334F) are shown in Fig|5| Maser emission was only 
detected in one source, 345.010+1.792, with thermal emis- 
sion detected towards Orion KL and NGC6334F. The ob- 
servation of NGC6334F was at an offset position and is 
of poor quality compared to previously published spectra 
jCragg et al. 2001) , so we have not included it in Fig [5] 
Marginal narrow emission was detected near the velocity 
of the 107.0-GHz maser peak towards two other sources 
(323.740-0.263 & 339.884-1.259) each with a significance of 
approximately 4-cr. 



4 DISCUSSION 

As well as assessing the prevalence of masers in the 85.5-GHz 
6_2-7_i E and 86.6-GHz 72-63 A - transitions of methanol, 
this work sets constraints on maser modelling. According to 
models of maser pumping, the detection of methanol maser 
emission in a star formation region signals the presence of 
special physical conditions, probably associated with a par- 
ticular evolutionary phase. The common 6.6- and 12.1-GHz 
methanol masers are excited for a wide range of model pa- 
rameters and so do not tightly constrain the physical condi- 
tions within the region. In contrast, modelling suggests that 
many of the other, rarer methanol maser transitions are in- 
verted over a more restricted range of parameters. The 85.5- 
and 86.6-GHz transitions each predominates over a differ- 
ent range of physical conditions (Crag g~et al. 200"T| , making 
them potentially useful probes of gas temperature and den- 
sity. 

Table [3] summarises the flux densities in 7 methanol 
maser transitions of the 25 star forming regions cur- 
rently known to harbour 107.0-GHz maser emission 
llCaswell et al. 20001 fVal'tts et al. 19951 . For the 6.6-, 12. 1-, 
107. 0-, 108.8- and 156.6-GHz transitions the flux densities 
cited are at the velocity of the 107.0-GHz peak and have 
been determined by eye from published spectra. For many 
of the sources detected at 108.8 and 156.6 GHz, thermal 
emission with no apparent maser component covers the ve- 
locity range of the 107.0-GHz peak. Where this occurs the 
upper limit listed in the table is marked with an asterisk and 
has been taken to be half the amplitude of the thermal emis- 
sion at that velocity. This approach has also been used in 
cases where thermal emission was detected from other tran- 
sitions. The most noticeable feature of Table [3] is that the 
85.5-, 86.6-, 108.8- and 156.6-GHz transitions are all uncom- 
mon (in the presence of 6.6-, 12.1- and 107.0-GHz masers). 
By investigating the range of physical conditions favouring 
each of these weaker transitions, we are able to distinguish 
between common and uncommon conditions within high- 
mass star formation regions. 



4.1 Variability in 345.010+1.792 

The 85.5- and 86.6-GHz methanol masers in 345.010+1.792 
were discovered by Cragg et al. ( 20(13 • Our survey con- 
firms the exceptional nature of this source, which remains 
the prime example of masers at these frequencies. Compar- 
ison of spectra for the same transitions obtained at differ- 
ent epochs or with different telescopes allows us to examine 
variability in these transitions, although this is complicated 
by uncertainties in the velocity and flux density calibration. 
FigureEK shows the 85.5-GHz emission in 345.010+1.792 at 
3 different epochs covering a period of 6 years. The single 
peaked line profile appears to have narrowed slightly over 
the period, but it is not possible to quantify the variability 
as some is likely to be due to calibration uncertainties. Fig- 
ure Oz> shows the 86.6-GHz emission in 345.010+1.792 at 3 
epochs covering essentially the same time range as the 85.5- 
GHz observations. The 86.6-GHz transition shows a clear 
double peaked profile and comparison of the relative inten- 
sities of these two features shows significant changes over a 
6 year interval. The amplitude calibration of the 1996 obser- 
vations from Mopra was uncertain by as much as a factor of 
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Table 1. A list of sources observed at 85.5-GHz. For sources where emission was detected the values listed are those of Gaussian profiles 
fitted to the spectra. For sources where no emission was detected the value listed is 3 x the RMS noise level in the Hanning smoothed 
spectrum. Upper limits indicated with a t were observed offset by 20" in right ascension and have been multiplied by a further factor of 
3 to account for reduced sensitivity. 



Name 


R.A.(J2000) 


Dec.(J2000) 


Integration 


Flux Density 


Velocity 


FWHM 




(h m s) 


(° ' ") 


Time (min) 


(Jy) 


(km s- 1 ) 


(km s- 1 ) 


Orion KT, 


05-35-14 5 


-05-22-30 


57 


25.0 


9.2 


2.6 










11.8 


8.3 


9.4 


188 9464-0 886 


06 08-53 3 


+21-38-29 


38 


<2.1 






1 Q2 600—0 048 


06-1 2-54 


4-1 7-59-24 


57 


<1.8 








±0 . tJ J. . tJO . tJ 


—61:15:42 


76 








e OAS— 1 Qfi 


1 c vOO r »^ 4 

J. o . \j yj . <j 'J . ti 


iJO. iJO. iJ'J 


38 


<2.0 






323 740—0 263 


1 5'31 -45 5 


-56-30-50 


38 


<1.8 






327.120+0.51f 


15:47:32.7 


-53:52:38 


76 


<1.5 






328.808+0.633 


15:55:48.5 


-52:43:07 


76 


1.4 


-43.7 


2.2 


336.0f8-0.827 


16:35:09.3 


-48:46:47 


76 


<1.3 






339.884-1.259 


16:52:04.7 


-46:08:34 


114 


<1.1 






340.785-0.096 


16:50:14.8 


-44:42:26 


38 


<1.8 






345.003-0.223 


17:05:10.9 


-41:29:06 


38 


<2.0 






345.010+1.792 


16:56:47.6 


-40:14:26 


37 


9.5 


-22.1 


0.7 


345.504+0.348 


17:04:22.9 


-40:44:22 


38 


< 4.7t 






348.703-1.043 


17:20:04.1 


-38:58:31 


38 


<1.5 






NGC6334F 


17:20:53.4 


-35:47:01 


37 


6.5 


-8.3 


5.9 


353.410-0.360 


17:30:26.2 


-34:41:46 


38 


<2.1 






9.621+0.196 


18:06:14.7 


-20:31:32 


38 


<1.5 






12.909-0.260 


18:14:39.5 


-17:52:00 


34 


< 5.0 1 " 






23.010-0.411 


18 34 40.3 


-09:00:38 


38 


< 5. 1 1 " 






23.440-0.182 


18:34:39.2 


-08:31:24 


19 


7.5 


100.7 


12.1 


35.201-1.736 


19:01:45.5 


+01:13:35 


91 


<1.1 







Table 2. A list of sources observed at 86.6-GHz. For sources where emission was detected the flux density and velocity are the values 
at the peak. For sources where no emission was detected the value listed is 3 x the RMS noise level in the Hanning smoothed spectrum. 
Upper limits indicated with a t were observed offset by 20" in right ascension and have been multiplied by a further factor of 3 to account 
for reduced sensitivity. 



Name 


R.A.(J2000) 


Dec.(J2000) 


Integration 


Flux Density 


Velocity 


FWHM 




(h m s) 


(° ' ") 


Time (min) 


(Jy) 


(km s _1 ) 


(km s _1 ) 


Orion KL 


05:35:14.5 


-05:22:30 


19 


15.8 


8.8 


2.1 










12.3 


8.1 


6.1 


188.946+0.886 


06:08:53.3 


+21:38:29 


34 


<2.4 






192.600-0.048 


06:12:54.0 


+17:59:24 


34 


<2.2 






310.144+0.760 


13:51:58.5 


-61:15:42 


38 


<2.1 






318.948-0.196 


15:00:55.4 


-58:58:53 


38 


<2.2 






323.740-0.263 


15:31:45.5 


-56:30:50 


76 


<1.9 






327.120+0.511 


15:47:32.7 


-53:52:38 


38 


<2.2 






339.884-1.259 


16:52:04.7 


-46:08:34 


34 


<2.4 






340.054-0.244 


16:48:13.9 


-45:21:44 


34 


<2.3 






340.785-0.096 


16:50:14.8 


-44:42:26 


34 


<1.9 






345.010+1.792 


16:56:47.6 


-40:14:26 


28 


16.4 


-22.0 


0.3 










10.1 


-21.2 


0.9 


345.504+0.348 


17:04:22.9 


-40:44:22 


31 


<2.4 






348.703-1.043 


17:20:04.1 


-38:58:31 


38 


<2.7 






NGC6334F 


17:20:53.4 


-35:47:01 


38 


4.0+ 






353.410-0.360 


17:30:26.2 


-34:41:46 


38 


<2.0 






9.621+0.196 


18:06:14.7 


-20:31:32 


38 


< 6.6t 






12.909-0.260 


18:14:39.5 


-17:52:00 


38 


< 8.6+ 






23.440-0.182 


18:34:39.2 


-08:31:24 


38 


<4.1 







Orion KL 
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328.808+0.633 





-50 -40 -30 

velocity w.r.t. LSR (kms+ 

NGC6334F 




velocity w.r.t. LSR (kms x ) 



velocity w.r.t. LSR (krns 1 ) 



23.440-0.182 




Figure 1. Spectra of the sources detected at 85.5 GHz. With the exception of 345.010+1.792 all spectra have been Hanning smoothed. 



2 jCragg et al. 200l| and comparison of the 86.6-GHz spec- 
trum of Orion KL from 1996 with that from the current 
observations shows the 1996 peak flux density to be lower 
by a factor of 1.8. The emission from Orion is thought to be 
quasi-thermal in origin and should therefore not vary, which 
suggests that the flux densities from the 1996 Mopra 86.6- 
GHz observations should be scaled in amplitude by a factor 
of 1.8. Note that all the spectra in Fig. [3]have been aligned 
with the peak emission velocity in observations made in Oc- 
tober 2002 using the newly upgraded Australia Telescope 



Compact Array (ATCA) at 3mm (Minier, Burton & Wong, 
in prep) . This required small, but significant shifts in some of 
the spectra, as detailed in the Figure caption, which can be 
attributed to an error in the AOS frequency calibration for 
the SEST 2000 observations, and intermittent errors in the 
velocity scale which are known to effect Mopra observations. 

The very narrow line width, offset velocity from the 
thermal emission in other transitions towards this source 
and clear temporal variability provides conclusive evidence 
that the 86.6-GHz emission in 345.010+1.792 is a maser. It 
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Orion KL 



345.010+1.792 




10 20 

velocity w.r.t. LSE (kms -1 ) 




-30 -20 -10 

velocity w.r.t. LSR (kms -1 ) 



Figure 2. Spectra of the sources detected at 86.6 GHz 



Table 3. The flux density at the velocity of the peak in the 107.0-GHz transition of the 6.6-, 12.1- , 85.5-, 86.6-, 107.0-, 108.8- & 156.6-GHz 
methanol transitions towards all star formation regions with known 107.0-GHz methanol masers. The 85.5- and 86.6-GHz data is taken 
from this work, except where indicated by a * . The information for the other transitions has been taken from the literature. Where upper 
limits are quoted they are 3 times the RMS noise level in the spectra. Transitions where thermal emission is detected are indicated with 
a * and the upper limited listed is 50% of the flux density of the thermal emission at the 107.0-GHz peak velocity. References : l=Batrla 
et al. <l987t : 2=Caswell et al. <l99^al ; 3=Caswe ll et al. 11995b! ; 4=Caswell et al. I2TO1 : 5=Cragg et al. lISUoTl ; 6=Koo et al. tT5ggl ; 
7=Mehringer, Zhou & Dickel ITW7V 8=Menten I1991bl : 9=Minier & Booth JOTSJ; 10=Slysh, Kalenskii & Val'tts (JTOSJ; ll=Sutton 
et al. |2T)0TI : 12= Val'tts et al. 11^951 : 13= Val'tts et al. llTTOl : 14=Walsh et al. tTWl 



Source 


Velocity 
(km s ) 


6.6-GHz 

(Jy) 


12.1-GHz 

(Jy) 


85.5-GHz 

(Jy) 


Flux Density 
86.6-GHz 107.0-GHz 

(Jy) (Jy) 


108.8-GHz 

(Jy) 


156.6-GHz 

(Jy) 


references 


W3(OH) 


-43.3 


3000 


600 


< 0.7 1 " 


6.7+ 


72 


<0.6 


<9 


8,10,11 


188.946+0.886 


10.9 


495 


235 


<2.1 


<0.8 


15.5 


<4.8 


<2 


2,3,4,5,13 


192.600-0.048 


4.2 


72 


<0.4 


<1.8 


<2.2 


5.8 


<4.5 


<3 


2,3,4,13 


310.144+0.760 


-56 


130 


114 


<1.4 


<2.1 


23 




<3 


4 


318.948-0.196 
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should also be noted that the strongest emission from the 
85.5- and 86.6-GHz transitions are coincident in velocity to 
within 0.1 km s _1 . The recent imaging of the 85.5- and 86.9- 
GHz methanol masers in 345.010+1.792 by Minier et al. 
has shown them to be spatially coincident at the arcsec- 



ond level at the same position (aj2ooo=16:56:47.56, <5j20oo=- 
40:14:26.3) as the 6.6-, 12.1- and 107.0 GHz masers. 
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Figure 3. (a) Spectra of the 85.5-GHz maser emission in 
345.010+1.792 at 3 different epochs. The velocity of the spec- 
trum taken at SEST in 2000 has been shifted by -0.3 km s _1 
and that of the spectrum from Mopra in 2002 has been shifted 
by -1.3 km s^ 1 . (b) Spectra of the 86.6-GHz maser emission in 
345.010+1.792 at 3 different epochs. The velocity of the spec- 
trum taken at Mopra in 1996 has been shifted by 0.3 km s _1 and 
that of the spectrum from SEST in 2000 has been shifted by -0.3 
km s^ 1 . The dotted line at -22.1 km s — 1 shows the velocity of 
the peak of the 85.5-GHz emission. 

4.2 Modelling 

The observations reported here, combined with the search of 
Minier & Booth (2002), demonstrate that class II methanol 
masers at 85.5 and 86.6 GHz are weak (< 14 Jy), and rarely 
detected at the sensitivity of the current experiment. We 
can use the results of maser modelling calculations to distin- 
guish between conditions in those few sources which display 
masers at these frequencies, the remaining target sources 
with methanol masers at both 6.6 and 107.0 GHz, and the 
majority of 6.6-GHz methanol maser sources which lack mil- 
limetre wavelength maser counterparts. 

We have used the model of Sobolev & Deguchi 
119941 to map out the prevalence of the various class II 
methanol masers as a function of the model parame- 
ters representing physical conditions. Calculations involv- 



ing this model were previously reported by Sobolev et al. 
jlQflfral H997pV Cragg et al. EHHTl, Su tton et al. j2~Uult . 
and Cragg, Sobolev & Godfrey (l2U02t . Altogether 1131 
energy levels of methanol (with torsional states v t = 
0, 1, 2 for both A and E symmetry species) were included 
jMekhtiev, Godfrey fc Hougen 1999} . The maser pumping 
is governed by infrared radiation from warm dust at temper- 
ature Td, which promotes the methanol molecules to their 
torsionally excited states. We assumed a dust filling factor 
of 0.5 and opacity (u/10 13 ) 2 , where v is the frequency in 
Hz. This is consistent with the observations of De Buizer, 
Pifia & Telesco 1200011 that the optical depth at 9.7 fim 
of mid-infrared emission associated with methanol masers 
can be as large as 1. Radiative transfer was treated in the 
large velocity gradient (LVG) approximation, augmented by 
a beaming factor (T x — 10, defined as the ratio of opti- 
cal depths parallel and perpendicular to the line of sight. 
This represents the assumed elongation of the maser region 
towards the observer. In this treatment the maser amplifi- 
cation and hence the brightness is governed by e _1 A r M/AT/, 
representing the column density of methanol along the line 
of sight divided by the line width of the emission. In terms 
of observable parameters, multiplying Nm/AV by the typ- 
ical line width of a methanol maser (0.5 km s _1 ) and the 
beaming factor gives the methanol column density. The pa- 
rameter Nm/ AV is referred to as the specific column density 
of methanol. The masers develop in gas of temperature Tk 
and density nn, where the methanol molecules undergo col- 
lisions with molecular hydrogen. Collisional excitation rates 
are based on the model of Peng & Whiteoak (19931, which 
uses propensity rules derived from a small number of labo- 
ratory measurements. The masers amplify radiation from a 
background ucHu region, here assumed to have brightness 
temperature 10 4 (1 - exp[-(1.2 x 10 10 /V) 2 ]) K, and geomet- 
rical dilution 0.002. The model produces methanol masers 
in the absence of a background ucH II region, but with lower 
brightness temperature in the centimetre wavelength masers 
than when the ucHu region radiation is included. 

Figure 2] displays the results of model calculations in 
which Td, Tk, tih and Nm/AV are systematically varied. 
Each plot shows the result of varying two of these param- 
eters, while the other two are held constant, representing 
selected planes through the parameter space. Figure 0] dis- 
plays a single contour for each of the maser transitions listed 
in Table El (6.6, 12.1, 85.5, 86.6, 107.0, 108.8, 156.6 GHz). 
The contour value is set at brightness temperature 10 6 K. 
Thus the area enclosed by each contour represents combina- 
tions of model parameters which give rise to masers at that 
frequency which are strong enough to be readily observed. 
Note that the exponential maser amplification leads to a 
very rapid rise in maser brightness, with the 6.6-GHz maser 
attaining brightness temperatures > 10 11 K. This style of 
plot is useful for identifying regimes where different masers 
switch on simultaneously. Such plots can also be used to de- 
fine limits on parameter values, while recognising that they 
represent only a partial exploration of the many dimensional 
phase space of the model. 

The top panel of Fig. [I] shows the behaviour of the 
model as dust and gas temperatures are varied, with the gas 
density fixed at riH = 10 6 ' 5 cm -3 and methanol specific col- 
umn density fixed at Nm/AV = 10 12 cm _3 s. Masers appear 
in the upper part of the plot, but disappear if the dust tern- 
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perature is too low (insufficient pumping). All the masers 
investigated here, with the exception of the 86.6-GHz line, 
are more readily excited at lower gas temperatures. Under 
the conditions shown, the 6.6- and 12.1-GHz masers appear 
together, requiring T d > 60 K. The 107.0-, 156.6- and 108.8- 
GHz masers require Td > 125 K, and progressively greater 
dust temperatures in warmer gas. The 85.5-GHz maser is 
confined to low gas temperature, while the 86.6-GHz maser 
is favoured by higher gas temperatures, with the two appear- 
ing together only for dust temperatures Td > 200 K and gas 
temperatures Tk < 40 K. De Buizer et al. (,20003 measure 
dust temperatures in the range of 100-200K for the mid- 
infrared sources associated with strong 6.6-GHz methanol 
masers, consistent with the general findings of the modelling. 

The lower panels of Fig. [I] show the behaviour of the 
model as methanol column density and gas density are var- 
ied, with the dust temperature fixed at Td = 150 K, for gas 
temperatures fixed at Tk = 30 and 130 K in the left and 
right panels respectively. Masers appear above the thresh- 
old in the upper left parts of the plots, but disappear if the 
methanol column density is too low (insufficient amplifica- 
tion) or the gas density is too high (thermalisation of popu- 
lations). The 6.6-GHz masers appear over the widest range 
of conditions, and are accompanied by masers at 12.1 GHz 
except at the lowest methanol column densities. The 107.0- 
, 156.6- and 108.8-GHz masers turn on in sequence as the 
methanol column density is increased, but do not extend to 
such high gas density as the centimetre wavelength masers. 
The 85.5-GHz maser is confined to low gas temperature and 
intermediate gas densities, while the 86.6-GHz maser has 
a somewhat complementary distribution, with the two ap- 
pearing together only for large methanol column density. 

The majority of class II methanol maser sources are 
characterised by centimetre wavelength maser emission at 
6.6 and/or 12.1 GHz, without maser counterparts at mil- 
limetre wavelengths. Points marked 'A' in Fig. 2] are exam- 
ples of model conditions producing masers above the thresh- 
old at 6.6 GHz only, while points marked 'B' produce 6.6- 
and 12.1-GHz masers simultaneously. The calculations sug- 
gest that these sources have dust temperatures 60 — 150 K, 
but do not place limits on their gas temperatures. The gas 
density is < 10 s cm -3 , with methanol specific column den- 
sity > 10 10,5 cm~ 3 s. It is likely that the ~ 50 percent of 
these sources which are detectable at 6.6 GHz but not at 
12.1 GHz have lower values of methanol column density than 
the sources exhibiting both masers. The few sites known to 
exhibit masers at 12.1 GHz only are probably at the high 
gas density limit of the maser range. 

In the following we assume that when masers at differ- 
ent frequencies are detected at matching velocities, they are 
simultaneously excited. Studies at high spatial resolution are 
required to confirm this by establishing whether or not the 
different masers are positionally coincident. In comparison 
with the centimetre wavelength maser sources, the 25 target 
sources which also have detectable masers at 107.0 GHz are 
likely to have slightly warmer dust temperatures, or to have 
slightly greater methanol column density, but do not exceed 
gas density 10 7 ' 5 cm -3 . Points marked 'C in Fig. [I]are ex- 
amples of conditions producing masers above the threshold 
at 107.0 GHz, but not at the other millimetre wavelengths, 
as was observed for 18 of the 25 sources in Table [3] The 
3 sources with masers at 156.6 GHz (but not at 85.5, 86.6 



or 108.8 GHz) continue the trend to greater dust temper- 
ature and/or methanol column density. Points marked 'D' 
in Fig.0]are examples of this combination. The remaining 4 
sources with 85.5-, 86.6- and/or 108.8-GHz masers represent 
outliers in the population of maser sites. 

We can estimate gas temperatures for these latter 4 
sources: first roughly from the regions in Fig. 2] with ap- 
propriate combinations of masers, and subsequently with 
a more quantitative approach to fitting model calculations 
to observations. For example, W3(OH) with masers at 6.6, 
12.1, 107.0 and 86.6 GHz (but not 85.5, 108.8 or 156.6 GHz) 
is likely to be an example of a warm gas source, with con- 
ditions resembling those of points 'E' in Fig. 0] This is con- 
firmed by detailed modelling by Sutton et al. H2001H which 
obtained T k = 150 K. The exceptional source 345.010+1.792 
shows maser emission in all 7 transitions, and if we assume 
these to coincide then Fig. [I] suggests a cool gas temper- 
ature, warm dust temperature and large methanol column 
density, as exemplified by the points marked 'F'. Detailed 
modelling by Cragg et al. 1200H confirms this (Model B 
of that paper), but predicts accompanying maser action at 
94.5 GHz which was not detected. 

The sources 328.808+0.633 and 9.621+0.196 exhibit 
masers at 6.6, 12.1, 107.0 and 85.5 GHz (but not at 86.6, 
108.8 or 156.6 GHz), although in both cases the 85.5-GHz 
maser is particularly weak. Fig. 2] suggests a cool gas tem- 
perature, with conditions resembling those of point 'G'. We 
have undertaken new fitting of the model to observations, 
along the lines described in Cragg et al. H2001H and Sutton 
et al. (2001}. Table 2] displays a comparison between the 
maser observations in 328.808+0.633 and 9.621+0.196 and 
examples of cool gas models (Tk = 20 K) which partially 
account for them. The models illustrated in Table 21 pre- 
dict strong masers at 6.6 and 12.1 GHz and weak masers 
at 85.5 and 107.0 GHz. The detection of very weak maser 
action at 85.5 GHz is consistent with these models, as are 
the nondetections at 86.6, 108.8 and 156.6 GHz. De Buizer 
et al. ( 2000 ) estimate the dust temperature in the region 
of 328.808+0.633 to be approximately 125K, in good agree- 
ment with the model in Table [I] However, they were unable 
to estimate a dust temperature for 9.621+0.196 due to the 
non-detection of mid-infrared emission towards this source. 
In both cases the calculated 12.1-GHz maser intensity is 
considerably greater than that observed. 

In order to compare observed flux densities with calcu- 
lated brightness temperatures, we must make an assumption 
about the size of the maser region. For the examples shown 
in Table |31 the size assumed is 60 milliarcseonds (mas) for 
328.808+0.633 and 20 mas for 9.621+0.196, chosen to match 
the flux density observed at 6.6 GHz to the brightness tem- 
perature calculated for these particular models. The detailed 
study of milliarcseconds structures of methanol masers by 
Minier et al. (I2002H found that at both 6.6 and 12.2 GHz the 
masers consist of a compact core (of the order of 0.5-3 mas in 
size) surrounded by a more extended halo with dimensions 
approximately an order of magnitude larger. The sizes we 
have assumed are therefore consistent with the dimensions 
of the halo component of the masers. One of the sources ex- 
amined in detail by Minier et al. {2002) was 9.621+0.196 at 
12.2 GHz. For the emission at -0.5 km s _1 they measured a 
halo size of > 50 mas and a core of 2 mas. 

The models presented above can account for the pres- 
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ence of different maser line combinations in the different 
sources, but do not agree quantitatively with all observa- 
tions, as Table 0] illustrates. This may be due to insufficient 
exploration of the parameter space, or to deficiencies in the 
model. While the maser pumping is dominated by radiative 
processes, the 85.5- and 86.6-GHz masers in particular are 
sensitive to the gas density and temperature, and therefore 
to the collisional excitation rates, which are not accurately 
known for methanol. Some calculations with nonselective 
collisional excitation rates by Sobolev et al. (If 997bB and 
Cragg et al. ((2001) illustrate how the maser intensities at 
high gas densities can be affected by changes to the colli- 
sion model. Until better collisional excitation rates become 
available for methanol, the density estimates from maser 
modelling will be subject to order-of-magnitude uncertain- 
ties. 

The examples above illustrate how the combination of 
methanol maser transitions observed in a particular source 
can be used to identify plausible excitation conditions. This 
approach is particularly valuable for sources which exhibit 
masers in several transitions simultaneously. The majority 
of class II methanol maser sources, however, have been de- 
tected in only the 6.6- and/or 12.1-GHz transitions. While 
the model can identify minimum conditions of dust temper- 
ature and methanol column density required to excite these 
masers, together with an upper limit to the gas density, ob- 
servations of these two lines alone are not enough to fully 
characterise the maser sites in these high-mass star form- 
ing regions. Nevertheless, the upper limits on nondetected 
lines provide a powerful constraint on the conditions in such 
sources, particularly when data is available for several lines 
sensitive to complementary factors. As further surveys are 
undertaken at a variety of methanol maser frequencies, it 
should be possible to refine quantitative estimates of physi- 
cal conditions, both in the exceptional sources with masers 
at many frequencies, and in the others which are more typ- 
ical of high-mass star formation regions. 



5 CONCLUSIONS 

We have undertaken a search for the 85.5- and 86.6-GHz 
transitions of methanol towards 22 southern sources which 
exhibit 1 07-GHz methanol masers. These observations show 
that maser emission from both of these transitions is rare. 
We have examined the implications of these findings in 
the context of the Sobolev & Deguchi model for methanol 
masers and find that our results are consistent with 

(i) The majority of class II methanol maser sources, which 
show methanol maser emission at only 6.6 and 12. 1 GHz, 
have dust temperature Ti = 60 — 150 K, gas density nn < 
10 s cm -3 , and methanol specific column density Nm/AV > 
10 105 cm" 3 s. 

(ii) The minority of sources which also show 107.0-GHz 
methanol maser emission have slightly higher Td or slightly 
larger Nm/AV than those which don't. 

(iii) The four known sources which also show 156.6-GHz 
methanol maser emission continue this trend to slightly 
larger T d and N M /AV. 

(iv) The four known sources with 85.5- and 86.6-GHz 
methanol maser emission represent outliers in the popu- 
lation of maser sites, with 85.5-GHz maser emission char- 



acteristic of low gas temperature (Tk = 20 K), 86.6-GHz 
maser emission characteristic of high gas temperature (Tk = 
150 K), and both present simultaneously only if Td > 200 K 
or Nm/AV > 10 125 cm~ 3 s. 

A comparison of the 85.5- and 86.6-GHz methanol 
maser spectra observed towards 345.01 0+1 .792 with pub- 
lished and unpublished spectra from earlier epochs shows 
clear variability, particularly for the 86.6-GHz transitions. In 
combination with the narrow line width of the emission and 
the offset of the 85.5- and 86.6-GHz emission from the ve- 
locity of the thermal emission for other molecular species in 
345.010+1.792 this presents conclusive evidence that these 
transitions are masing in this exceptional source. 
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the parameter values are Tjj = 20 K, = 
are as described in Section l4.2l of the text. 
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85.5- and 86. 6- GHz methanol maser emission 11 
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Figure 4. Contour diagrams showing regions where selected class II methanol masers become active, as a function of 2 model parameters. 
Contours are labelled with the maser frequency in GHz. For each maser a single contour representing brightness temperature 10 6 K 
is drawn, with no masers appearing above this threshold in the lower regions of the plots. Top panel shows behaviour as T^ and TJj 



are varied, when nji 



10" 



cm i and TVm/AV = 10 cm s s. Bottom panels show behaviour as njj and JVm/AV are varied, when 



Tj = 150 K for Tit = 30 K (left panel) and T^ = 130 K (right panel). Points marked 'A' to 'C illustrate different combinations of masers 
present simultaneously, as discussed in Section l4.2l of the text. 



